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Abstract 

Background and Aims: kzaty\\opr\ne (AZA) is widely used for the treatment of inflammatory bowel disease (IBD) patients. 
AZA is catabolized by thiopurine S-metliyltransferase (TPMT), wliicli exhibits genetic polymorphisms. It has also been 
reported that 5-aminosalicylic acid (5-ASA) inhibits TPMT activity, and that increased 6-thioguanine nucleotide (6-TGN, a 
metabolite of AZA) blood concentrations result in an increased number of ADRs. In this study, single nucleotide 
polymorphisms (SNPs) related to differential gene expression affecting AZA drug metabolism in combination therapy with 
5-ASA were examined. 

Methods:lo identify genetic biomarkers for the prediction of 6-TGN blood concentration, ExpressGenotyping analysis was 
used. ExpressGenotyping analysis is able to detect critical pharmacogenetic SNPs by analyzing drug-induced expression 
allelic imbalance (EAI) of premature RNA in HapMap lymphocytes. We collected blood samples on 38 patients with 
inflammatory bowel disease treated with AZA and corroboration of the obtained SNPs was attempted in clinical samples. 

Results: k large number of SNPs with AZA/5-ASA-induced EAI within the investigated HapMap lymphocytes was identified 
by ExpressGenotyping analysis. The respective SNPs were analyzed in IBD patients' blood samples. Among these SNPs, 
several that have not yet been described to be induced by AZA/5-ASA were found. SNPs within SLC38A9 showed a 
particular correlation with patients' 6-TGN blood concentrations. 

Conclusions: Based on these results, ExpressGenotyping analysis and genotyping of patients appears to be a useful way to 
identify inter-individual differences in drug responses and ADRs to AZA/5-ASA. This study provides helpful information on 
genetic biomarkers for optimized AZA/5-ASA treatment of IBD patients. 
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Introduction 

The immuiie-modulatiiig thiopurines, 6-mercaptopurine (6-MP) 
and its pro-drug azathioprine (AZA), are widely used in 
inflammatory bowel disease (IBD) treatment [1—3]. AZA is a 
pro-drug of 6-MP carrying a methyl imidazole group attached to a 
sulfur carbon of 6-MP and metabolized to 6-MP by cleavage 
require the action of gluthatione S transferase. For the mainte- 
nance of remission in IBD, AZA and 6-MP are the first-line 
immunomodulators [4,5]. Both AZA and 6-MP must be 
extensively metabolized before the pharmacologically active 



metabolites, 6-thiogTaaninenucleotides (6-TGN), are generated. 
The clinical efficacy of 6-MP has been reported to be correlated 
with erythrocyte levels of 6-TGN [6,7]. The usefulness of 
monitoring 6-TGN blood concentrations in patients receiving 
AZA and 6-MP has also been suggested [8-10]. Although these 
thiopurine derivatives are considered to be relatively safe for 
maintenance therapy, several studies have reported discontinua- 
tion of thiopurine derivatives in up to 50% of patients during long- 
term therapy, mainly due to the development of adverse drug 
reactions (ADRs) [4, 1 1 , 1 2] . AZA is catabolized by thiopurine S- 
methyltransferase (TPMT), which exhibits single nucleotide 
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polymorphisms (SNPs). These SNPs have been reported to 
contribute to bone marrow suppression as ADRs [9,13,14]. 
Patients with mutant alleles are reported to exhibit lower 

enzymatic activit}' of TPMT and to be more likely to show bone 
marrow suppression. In the Japanese population, TPMT A719G is 
the most common allele [15-18] and is associated with low TPMT 
enzyme activity [19,20]. However, in the Japanese population, the 
frequency of TPMT mutations is reported to be lower (2%) than 
that reported from Western countries [15-17,21,22]. It has also 
been reported that 5-aminosalicylic acid (5-ASA) inhibits TPMT 
activity, and that an increased 6-TGN blood concentration results 
in an increased number of ADRs [3]. However, the mechanism of 
ADRs induced by AZA has not yet been completely understood. 
Therefore, to identify new genetic bi()mark(^rs to predict 6-TGN 
blood concentrations before the administration of thiopurine is 
very valuable and important. 

In this study, SNPs related to differential gene expression 
affecting AZA drug metabolism in AZA/5-ASA combination 
therapy were investigated. This genetic information is important to 
predict the 6-TGN blood concentrations of IBD patients treated 
with AZA/5-ASA. Recendy, ExpressGenotyping analysis has been 
developed to explore the SNPs affecting gene expression changes 
of each allele (these gene expression changes are called Expression 
Allelic Imbalance (EAI)) after drug exposure (Figure 1). It is clear 
from a previous report that the concept of allelic variation such as 
this EAI is also important [23-25]. Such SNPs may be associated 
with inter-individual differences in drug response and ADRs. 
Therefore, ExpressGenotyping analysis can be used as a ntnx-l 
technology to detect critical genetic biomarkers that predict 6- 
TGN blood concentrations of IBD patients. 

Materials and Methods 

Patients 

Patients with established IBD receiving AZA were enrolled 
between April 1^''' 2009 and August 31*^^ 2009. All patients were 
Japanese persons attending the gastroenterology outpatient clinic 
at the Hospital of th(- Kyoto Prc-icctural University of Medicine. 
The characteristics of the enrolled patients are shown in Table 1 . 
Of the 38 patients, 35 were also treated with 5-ASA. Clinical 
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information related to drug administration for these patients is 
shown in Table 2. The protocol of this study was approved by the 
Ethics Committee of the Kyoto Prefectural University of Medicine 
and written informed consent was obtained from all patients 
before the original procedures were conducted. This study was 
performed in accordance with the ethical principles associated 
with the Declaration of Helsinki and registered in the University 
Hospital Medical Network Clinical Trials Registry (UMIN 
000010197). 

Measurement of 6-TGN 

The blood 6-TGN concentration was measured by Mitsubishi 
Chemical Medience Corporation using HPLC (InertsU ODS-3 
(3.0 x100 mm, 4 |j^m). This is the modified method from the 

previous report [23]. 

Exposure of HapMap lymphocytes to drugs 

AZA and 5-aminosalicylic acid (5-ASA) were obtained from MP 
Biomedicals, Inc. (Tokyo, Japan). They were dissolved in RPMI 
medium 1640 (Invitrogen, Van Allen Way Carlsbad, CA, USA) 
including 15% fetal bovine serum (Invitrogen) as vehicle for in 
vitro experiments. Suspensions of 30 different HapMap lympho- 
cytes from Japanese subjects (Coriell, Camden, INf, USA), at a 
concentration of 2.5x10® cells in 0.5 mL of RPMI medium 1640 
(Invitrogen) including 15% fetal bovine serum, were incubated 
with 4.5 mL of AZA/5-ASA (100 (xM AZA and 10 mM 5-ASA). 
The suspensions were sampled after 24 hours and used for the 
following assay. As controls, the respective HapMap lymphocytes 
were treated in the same way, but only with the vehicle RPMI 
medium 1640 including 15% fetal bovine serum (no drugs). 

mRNA expression analysis 

Total RNA was extracted from the AZA/5-ASA-treated and 
control cells using TRIzol Reagent (Invitrogen). Total RNA 
concentration was determined using a NanoDrop ND-1000 
instrument (NanoDrop Technologies, Wilmington, DE, USA), 
and RNA quality was checked with an Agilent 2100 Bioanalyzer 
(Agilent Technologies, Santa Clara, CA, USA). Starting with 
250 ng of total RNA, GeneChip U133 Plus 2.0 Array analysis 
(Asymetrix, Santa Clara, CA, USA) was conducted following the 
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Figure 1. ExpressGenotyping analysis. ExpressGenotyping analysis can detect critical pharmacogenetic SNPs by analyzing drug-induced 
expression allelic imbalance (EAI) of premature RNA in HapMap lymphocytes. This new analysis technique has the potential to accelerate the 
realization of personalized medicine, which is closely associated with inter-individual differences in drug response and ADRs. 
doi:1 0.1 371/journal.pone.0095080.g001 
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Table 1. Characteristics of enrolled patients. 




Total (n) 


38 


Sex (n; male/female} 


29/9 


Age (years; median±range) 


44.0±10.6 


Disease (n; UC/CD/BD) 


21/14/3 


Location of UC patients (n = 21) 


Pancolitis 


19 


Distal colitis 


2 


Location of CD patients (n = 14) 


Colon 


4 


Colon and ileum 


6 


Ileum 


4 


Duration of IBD (years; median±range) 


12.8 ±8.97 


UC: ulcerative colitis, CD: Crohn's disease, BD: Behcet's disease. 




doi:1 0.1 371 /journa!.pone.0095080.t001 





manufacturer's protocol, using the GeneChip 3'IVT Express Kit 
(Affymetrix) and the GeneChip Hybridization Wash and Stain Kit 
(AflFymetrix). The raw expression signal obtained through a 
GeneChip Instrument System (Affymetrix) was normalized by 
scaling of the target signal to 100. For further data processing, 
probes were selected showing a P call in the detection step of the 
U133 Plus 2.0 Array among these selected probes; any probe set 
showing a &2-fold increase or decrease (AZA/5-ASA- treated vs. 
vehicle) of gene expression for a given sample was considered a 
significant change. 

SNP 6.0 array analysis 

DNA was extracted from the 30 different untreated HapMap 
lymphocytes using a QIAamp DNA mini Kit (QIAGEN), and the 
quality and concentration of the DNA samples were determined 
with a NanoDrop ND-1000 (NanoDrop Technologies). On the 
other hand, DNA of IBD patients was extracted with FALCO 
Biosystems (Kyoto, Japan). Starting with 500 ng of DNA, the 6.0 
SNP Mapping Array assay was performed using the GeneChip 6.0 
SNP Mapping Assay Kit (AflFymetrix) following the manufacturer's 
protocol. Data for SNP typing were determined using the 
Genotyping Analysis Software (GTYPE) Version 4.1 of the 
GeneChip Instrument System (Affymetrix). Based on the 6.0 
SNP Mapping Array genotyping information, an AA, AB, or BB 
call was obtained. In this analysis, a Birdseed vl or v2 algorithm 
was included in the calculation [26]. 

MassARRAY analysis 

DNA of IBD patients was extracted with FALCO Biosystems. 
Genotyping of SNPs (rsl 142345, rsl800584, rs56161402, 
rs6921269, rs2842951, rs2842934, rsl 800460, rs25 18463, 
rs72552739, rsl800462, and rs3898137) in the TPMT gene was 
performed using the Sequenom MassARRAY technology platform 
(Sequenom, San Diego, CA, USA), with the iPLEX Gold Reagent 
Kit (Sequenom), SpectroCHIP ARRAYS, Clean Resin Kit 
(Sequenom), and HotStarTaq DNA Polymerase (QIAGEN, 
HUden, Germany) based on the single-base extension reaction. 
Locus-specific PGR primers and allele-specific detection primers 
were designed using the MassARRAY Assay Design 3.1 software 
(Sequenom). Allele detection was performed using a MassARRAY 
Analyzer Compact MALDI-TOF mass spectrometer (Sequenom), 



a MassARRAY Nanodispenser (Sequenom). 
RAY TYPER 3.4 software (Sequenom). 



and the MassAR- 



Correlation analysis of SNPs and 6-TGN risk ratios in IBD 
patients 

We difmed the risk ratio as the proportion of blood 6-TGN 
concentration and the dose of AZA intake (6-TGN concentration/ 
AZA dose). The correlations between the AA, AB, and BB 
genotypes and 6-TGN risk ratios of clinical data in IBD patients 
were analyzed. For the correlation analysis, the genotypes of each 
sample from 32 individuals were plotted on the X-axis as AA = 0, 
AB = 1, and BB = 2, and 6-TGN risk ratios for the corresponding 
genes were plotted on the Y-axis. These for which all AA, AB, and 
BB genotypes were observed in at least 2 samples were selected for 
subsequent analysis, and r^ was calculated for the combination of 
each SNP and 6-TGN risk ratio. With this correlation analysis, 
SNPs with r^>0.4 (p-values<0.01) with 6-TGN risk ratios were 
identified, and these were considered genetic biomarkers. 

ExpressGenotyping analysis 

Premature RNA expression analysis was performed using the 
ExpressGenotyping reactor and ExpressGenotyping analyzer 
(Japanese Patent Number: 4111985, Method for identifying a 
gene with allelic variation in gene expression (FILING DATA: 
December 15, 2005)). This new analysis technique is performed 
using the GIM algorithm [27-29]. The ExpressGenotyping 
reactor consists of a group of reaction steps where cDNA is 
synthesized preferentially from the premature RNA in the total 
RNA. The subsequent assay using the synthesized cDNA is 
essentially similar to the Affymetrix 6.0 SNP Mapping Array's 
assay, but instead of ordinary DNA typing, quantitative allele- 
specific expression information is obtained. The ExpressGenotyp- 
ing analyzer combines the information from the DNA and cDNA 
typing and detects variations of premature RNA expression from 
each allele at particular SNPs. Consequently, comprehensive 
premature RNA expression assays using the ExpressGenotyping 
reactor and ExpressGenotyping analyzer can detect genes with 
diflferential allelic expressions. 
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Table 2. Clinical Information of drug administration for patients. 









Dose of azathioprine 


Dose of 5-ASA 


Conc6ntr3tion of 6-TGN 




Sample number 


Catalog ID 


(mg/day) 


(mg/day) 


(pmol/8x10' RBCs) 


Risk ratios (6-TGN/AZA) 


1 


5842535 


75 


3600 


195 


2.60 


2 


5779006 


50 


2250 


474 


9.48 


3 


770916 


25 


1500 


423 


16.92 


4 


4872959 


25 


3000 


130 


5.20 


5 


5459355 


50 


3000 


287 


5.74 


6 


6232762 


50 


2250 


243 


4.86 


7 


2639084 


50 


2250 


164 


3.28 


8 


0156158 


50 


1500 


144 


2.88 


9 


5571916 


50 


3000 


50 


1.00 


10 


3281892 


25 


2250 


464 


18.56 


11 


2852065 


100 


3000 


353 


3.53 


12 


162980 


50 


2250 


572 


11.44 


13 


0777935 


100 


0 


253 


2.53 


14 


5445356 


100 


2250 


178 


1.78 


15 


6007418 


50 


3600 


152 


3.04 


16 


5294675 


25 


2250 


144 


5.76 


17 


4549910 


100 


3000 


87 


0.87 


18 


5999831 


25 


2400 


163 


6.52 


19 


2002116 


25 


2400 


54 


2.16 


20 


6158273 


50 


4000 


365 


7.30 


21 


6176417 


50 


3600 


293 


5.86 


22 


4255178 


50 


3600 


557 


11.14 


23 


6230023 


50 


3600 


99 


1.98 


24 


0185310 


75 


2000 


595 


7.93 


25 


0437466 


50 


3000 


488 


9.76 


26 


1830237 


50 


2250 


359 


7.18 


27 


1996388 


50 


3000 


146 


2.92 


28 


2240525 


50 


1500 


301 


6.02 


29 


2366881 


50 


3000 


287 


5.74 


30 


4782476 


50 


2000 


193 


3.86 


31 


4916624 


50 


3000 


201 


4.02 


32 


5478242 


100 


2250 


380 


3.80 


33 


5480939 


25 


3000 


292 


11.68 


34 


5615377 


100 


4000 


113 


1.13 


35 


5796199 


100 


0 


50 


0.50 


36 


5892802 


75 


3000 


366 


4.88 


37 


5906785 


50 


2000 


561 


11.22 


38 


6176174 


50 


0 


288 


5.76 



Clinical information related to drug administration for patients. 

Since the 5-ASA dose was less than 2000 mg/day in sample numbers 3, 8, 13, 28, 35, and 38, they were not used for data analysis in this study. 
doi:l 0.1 371 /journal.pone.0095080.t002 



Detection of drug-induced EAI by ExpressGenotyping 
analysis 

To identify genetic biomarkers involved in tlie AZA/5-ASA- 
response, ExpressGenotyping analysis was performed as follows: 1) 
The ExpressGenotyping data outputs of drug-treated and vehicle- 
treated cells were obtained as a text file; 2) The respective vehicle 
and drugs files were combined with probe ID; 3) Probes with a 
heterozygous (AB) call within a gene region (3UTR, 5UTR, CDS, 
exon, and intron) were selected; 4) The expression ratios of 



premature RNA (vehicle sample (A aUele/B allele or B allele/A 
allele) "vehicle EAI (max)", drug sample (A allele/B allele or B 
allele/A allele) "drug EAI (max)" were calculated; 5) SNPs with an 
EAI ratio [calculated as drug EAI (max)/vehicle EAI (max)] of S2 
were considered to show an EAI specific for drug exposure; and 6) 
EAI data of 30 test samples that were cell lines were collected. 
Next, AB heterozygous genotype analysis based on statistics was 
performed to identify responsive SNPs from the cells exhibiting the 
AB heterozygous genotype. SNPs with 3 or more heterozygous 
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Figure 2. Frequency distribution of 6-TGN levels with the treatment of AZA and 5-ASA. (A) 6-TGN levels of patients taking less than 
3000 mg/day of 5-ASA. (B) 6-TGN levels of patients taking more than 3000 mg/day of 5-ASA. Relationships between 6-TGN levels and MCV (C) and 
between 6-TGN levels and WBC (D). 
doi:l 0.1 371/journal.pone.0095080.g002 



calls among the 30 test samples were used for this analysis. In this 
method, SNPs showing p-values<0.01 in either analysis of covariance 
of expre.ssion values or the <-test of the difference of inter-allelic 
expression ratios in cells under stimulation with test compounds 
and also exhibiting inter-allelic expression ratios of s2-fold were 
added to the list of test compound-responsive SNPs. 

Comparative analysis between clinical and 
ExpressGenotyping data 

The genetic biomarkers analyzed in IBD patients' blood 
samples were correlated with drug-induced EAI obtained by 

ExpressGenotyping analysis. In this analysis, a coincidence 
between genetic biomarkers of clinical data (r^ = 0.4 or more, p- 
value = 0.01 or less) of IBD patients and drug-induced EAI values 
of S2 was detected based on the statistics of the AB heterozygous 
genotype analysis. 

Statistics 

Risk ratio was assessed for each genotype group of patients. The 
data were compared using correlation analysis and t test. Prizm 
version 5.0 d (GraphPad Software Inc, La Jolla, CA) and SciPy 
(statistical library of Python) (http:/ /www.scipy.org/) were used 
for all statistical analysis. All statistical significance was defined as 
P<.05. 



Results 

6-TGN blood concentration and 5-ASA dose 

The relationship between 6-TGN blood concentration and the 
5-ASA dose was analyzed. There was no positive correlation 
between 6-TGN blood concentration and low-dose 5-ASA 
(Figure 2a) or high-dose 5-ASA (Figure 2b) administration. 

As shown in Figure 2c, mean corpuscular volume (MCV) was 
highly correlated with 6-TGN blood concentration. The white 
blood cell (WBC) level showed a significandy negative correlation 
with 6-TGN levels (Figure 2d). 

mRNA expression analysis 

Gene expression changes of mRNA were observed in mRNA 
expression analysis of HapMap lymphocytes from Japanese 
subjects. The scatter plots (Figure SI) showed significant changes 
in mRNA expression for the 100 |xM AZA and 10 mM 5-ASA- 
treated vs. vehicle-treated HapMap lymphocytes for 24 hours. As 
the probe numbers of gene expression changed, by observing an 
average of 1257 up (change ^2 x)-regulated-probes (796 genes) 
and 815 down (change SO.5 x)-regulated-probes (621 genes), 
100 nM AZA and 10 mM 5-ASA treatment caused a marked 
effect on mRNA expression. GOTERM analysis of the Database 
by Annotation Visualization and Integrated Discovery (DAVID) 
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Figure 3. SLC38A9 SNPs of IBD patients. SLC38A9 SNPs were determined by correlation analysis with genotypes and 6-TGN risk ratios in IBD 
patients. 

doi:1 0.1 371 /journal.pone.0095080.g003 



showed the up-related terms (apoptosis, programmed cell death, 
and death etc.) and the down-related-terms (M phase, cell cycle 
phase, and cell cycle etc.) (http://david.abcc.ncifcrf.gov/). Based 
on the above results, exposure with 100 |xM AZA and 10 mM 5- 
ASA for 24 hours was selected as the experimental condition for 
the subsequent ExpressGenotyping analysis of 30 different 
HapMap lymphocytes (Figure S2). 

SNPs analysis in IBD patients 

Correlation analysis of SNPs and 6-TGN risk ratios in IBD 
patients with SNP 6.0 Array analysis detected 1286 SNPs (516 
genes) with r>0.4 (p-values<0.01) (Figure S3). The GOTERJvI 
analysis of the database by DAVID identified cell adhesion, 
biological adhesion, and cell projection morphogenesis, etc. The 
top most significant SNPs were genotypes for rs3846502, 
rs3761769, rsl6884434, and rs6897117 in SLC38A9 (Figure 3). 
SLC38A9 has been reported as a transporter gene. On the other 
hand, no correlations with SNPs in TPMT gene were obtained. 
Therefore, in order to perform a more detailed analysis of the 
TPMT gene, by MassARRAY analysis, the genotypes (AA, AB, or 
BB) of various TPMT SNPs were identified. TPMT is an enzyme 
that metabolizes AZA to 6-TGN, and some SNPs cause 
differences in metabolic activity. Four SNPs of TPMT have been 
analyzed (rsl800462, TPMT*2 [30-32]; rsl800460, TPMT*3B 
[33]; rsl 142345 "A719G", TPMT*3C [34]; rsl800584, TPMT*4 
[35]) that affect 6-TGN blood concentration and, therefore, 
hematological toxicity. No difference in the distribution of SNPs 
was found, and genotype effects of the risk allele on hematological 
toxicity were not seen in all clinical samples measured (Figure S4). 



Other taq SNPs (rs56161402, rs6921269, rs2842951, rs2842934, 
rs2518463, rs72552739 and rs3898137) in TPMT showed no 
differences in these genotypes, as with the SNPs noted above. 
These SNPs of rs2842951, rs2518463, and rs3898137 observed in 
different distributions did not show significant correlations 
(rs2842951 (r = 0.08014, p-values = 0.66281), rs2518463 
(r = 0.10855, p-values = 0.55427), and rs3898137 (r = 0.04124, p- 
values = 0.82268)). From the above results, differences in the 6- 
TGN risk ratios were found to not be effected by TPMT SNPs of 
IBD patients. 

Detection of drug-induced EAI 

For the detection of AZA/5-ASA-induced EAIs, 8373 SNPs 
(3338 genes) with an liAI ratio S2 in ExpressGenotyping analysis 
(Figuer S5) were selected, and the numbers of AZA/5-ASA- 
induced EAIs for each HapMap lymphocyte are shown in 
Figure 4a. The dispersion diagram of the AZA/5-ASA-induced 
EAIs in the AB heterozygous genotype analysis for each 
chromosome is shown in Figure 4b. The most significant genes 
showing &3 fold-change of EAI were the SNPs of chromosome 
11: GALNTL4 (rs2129523), chromosome 14: NPAS3 
(rsl2433161), chromosome 20: PHACTR3 (rs6015561) and 
chromosome 20: MACROD2 (rs81 16068). In addition, the AB 
heterozygous genotype analysis of ExpressGenotyping data 
revealed reliable AZA/5-ASA-induced EAIs (64 SNPs (55 genes). 
The results are shown in (Figure S6). GOTERM analysis of the 
database with DAVID identified the regulation of apoptosis, 
regulation of programmed cell death, and regulation of cell death, 
etc. in the AB heterozygous genotype analysis. These EAI-related 
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Figure 4. Number of drug-induced EAIs for 30 different HapMap lymphocytes by ExpressGenotyping analysis. (A) For the drug- 
induced EAIs, SNPs with an EAI ratio >2 were selected. Dispersion diagram of drug-induced EAIs for each chromosome. (B) The drug-induced EAI of 
fold-change in the dispersion diagram detected with AB heterozygous genotype analysis by ExpressGenotyping. Details of the AB heterozygous 
genotype analysis are described below. SNPs with the AB heterozygous genotype in >3 of 30 HapMap lymphocytes were used to search for SNPs 
responsible for producing the drug-induced intracellular allelic imbalance. These analyses were performed using the SNPs for which the AB 
heterozygous genotype was observed in at least 3 samples. 
doi:1 0.1 371 /journal.pone.0095080.g004 



SNPs on ExpressGenotyping analysis may be important as genetic 
biomarkers affecting gene expression changes of each allele in 
predicting 6-TGN blood concentrations in patients taking AZA/5- 
ASA medication. 

Comparative data analysis by clinical data and 
ExpressGenotyping analysis 

An attempt was made to identufy genetic biomarkers by 
combining 6-TGN risk ratios of clinical data and AZA/5-ASA- 
induced EAI of ExpressGenotyping analysis. In this comparative 
data analysis, SNPs in SLC38A9 showed a correlation with the 6- 
TGN risk ratios and AZA/5-ASA-induced EAI. The SNP 
rs6897117 of SLC38A9 was significant in the clinical data, and 
ExpressGenotyping analysis showed an EAI ratio &2 in GM18975 
(Figure 5). The EAI ratio of SLC38A9 showed that the gene 
expression level of allele B was reduced, and it may increase the 6- 
TGN risk ratios. Therefore, these SNPs identified by Express- 
Genotyping analysis might serve as genetic biomarkers for 
prediction of 6-TGN blood concentrations in patients taking 
AZA/5-ASA medication. 

Discussion 

In the present study, new SNPs on SLC38A9 as genetic 
biomarkers to predict blood 6-TGN concentrations in patients 
taking AZ A were identified. As the usage of AZ A for IBD patients 
is now increasing, and the efficacy of AZA to maintain remission 



has been reported to be superb, it is very important that 
complications do not occur with AZA. In the analysis of TPMT 
SNPs of IBD patients in this study, the difference in the 6-TGN 
risk ratios was found to not have the risk allele related to 
hematological toxicity. Compared to the previous genetic 
biomarker, the new biomarker, SNPs of SLC38A9, was more 
sensitive. 

In this study, 6-TGN concentration was not correlated with the 
dose of AZA treatment. Most of the studies have also reported no 
relationship between AZA dose and 6-TGN concentration [36]. 
This means that weight-based dosing of AZA does not guarantee 
therapeutic metabolite concentrations, and the efficacy of AZA 
metabolism might determine 6-TGN concentrations. As the 
number of WBCs and 6-TGN concentration showed a reverse 
correlation in the present study, prediction of blood 6-TGN 
concentration before AZA treatment will contribute to prevent 
leukopenia induced by 6-TGN. 

Solute carriers (SLCs) are the largest group of transporters in 
the human genome [37,38]. Many of the SLCs have important 
functions in cellular processes, including roles as transporters for 
amino acids and neurotransmitter cycling, while others are passive 
transporters or exchangers [38-40]. The SLCs have been 
categorized into at least 46 different families with varied 
biochemical properties [41], and two more families have recendy 
been identified [42]. In mammals, the SEC proteins form four 
major phylogenetic groups, a-, fi-, y-, and 5-groups, with proteins 
in each group having a common evolutionary origin. The f5-group 
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Figure 5. Results of drug-induced EAI for SNP (rs6897117) of 
SLC38A9. The SNP of SLC38A9 shows a consistent relationship 
between 6-TGN risl< ratios and drug-induced EAI. 
doi:1 0.1 371 /journal.pone.0095080.g005 

includes the SLC32, SLC36, and SCL38 families, and it is the 
second largest phylogenetic cluster of amino acid transporters. 

The SLC38 family consists of 1 1 members. The orphans 
SLC38A7-1 1 were only recently identified, and their substrate of 
transport has thus far not been characterized [43], while the other 
members, SLC38A1-6, are functionally classified as sodium- 
coupled neutral amino acid transporters, also known as SNATs. 
SNATl (SLC38A1, [44,45]), SNAT2 (SLC38A2, [46^9]), and 
SNAT4 (SLC38A4, [50,51]) have been further classified into 
system A, while SNAT3 (SLC38A3) and SNAT5 (SLC38A5) 
belong to system N transporters [52-54]. The system N/A 
classification depends on the functional properties and patterns of 
substrate recognition. SNAT6 (SLC38A6) [55] and the other more 
recently discovered members of the family have not been classified 
according to the N/A system. AH SNAT genes, except SNAT4 
[50], are expressed in the brain, with SNATl [44,56] and SNAT2 
[57] expressed both in astrocytes and neurons, while SNAT3 [52] 
and SNAT5 [58] are located in glial cells. No exact cellular 
localization of the protein or substrate profile has been reported 
for the orphan members of the SLC38 family. L-glutamine is a 
favored substrate for most of the members throughout the SLC38 
family, and it is suggested that the transporters are part of the 
glutamine-glutamate cycle in the CNS [52,54]. In this study, the 
genotype of SLC38A9 was found to be related to the blood 6- 
TGN concentration. The direct effect of SLC38A9 SNPs on 6- 
TGN blood concentration is still unclear, and further study will be 
necessary to reveal the pathway from SLC38A9 SNPs to the 
metabolism of AZA. 

ExpressGenotyping analysis developed by HaploPharma Inc. is 
an innovative technology to detect critical pharmacogenetic SNPs 
and to provide genetic biomarkers to analyze them. ExpressGen- 
otyping analysis enables the discovery of genetic biomarkers with a 
much smaller number of samples than with conventional 
approaches, such as genome-wide association study (GWAS), 



because these conventional approaches require complex proce- 
dures that usually need samples of hundreds or thousands of cases 
and a long time period, as well as being quite costly. In comparison 
to these approaches, ExpressGenotyping analysis can analyze the 
expression of each allele in a sample separately and, thus, reveal 
differences in allelic gene expression within the sample. Therefore, 
ExpressGenotyping analysis can be used to perform practical 
genetic biomarker discovery in samples with smaller numbers, at 
low cost. The basic functional principle of ExpressGenotyping is as 
follows. 1) Cell lines (HapMap lymphocytes) are treated with the 
drug of interest. 2) RNA and DNA are prepared from drug-treated 
and untreated cell lines. 3) Combined SNP typing/ expression 
profiling on Affymetrix GeneChip microarrays is performed. 4) 
Analysis of expression profiles (drug-treated vs. vehicle-treated) 
and drug- induced EAI is then done. 5) Candidate genetic 
biomarkers for drug response prediction are identified. 6) 
Candidate genetic markers are validated with sample material 
from clinical studies. These technologies combined allow the 
genome-wide detection of variations in gene expression of 
individual alleles within one sample and eliminate the requirement 
for extensive inter-sample statistical analyses to identify a 
correlation between gene expression and genotype. In this report, 
the genetic biomarkers were identified by combining clinical data 
of IBD patients with AZA/5-ASA-induced EAI obtained by 
ExpressGenotyping analysis. Many SNPs in different genes were 
identified as candidate genetic biomarkers showing the relation- 
ship biologically (GOTERM analysis of DAVID). Of the genetic 
biomarkers related to the clinical 6-TGN blood concentration, 
SNPs in SLC38A9 were the most significant. This SLC38A9 gene 
might have an important role as a transporter for various drugs. 
Therefore, this SNP identified by ExpressGenotyping analysis 
might serve as a genetic biomarker for prediction of the 6-TGN 
blood concentration in patients on AZA/5-ASA. 

These results show that ExpressGenotyping analysis is a useful 
way to identify individual differences in the incidence of drug 
response and ADRs, and it is expected to be a new method to 
accelerate the achievement of personalized medicine. In the 
future, the EAI EG method can be applied to the analysis of ex 
ante evaluation of immunosuppressive agents in the treatment of 
IBD, such as AZA. 

Supporting Information 

Figure SI Scatter plot analysis. In order to determine the 
experimental conditions for ExpressGenotyping analysis, scatter 
plot analysis was performed with the U133 Plus 2.0 Array. Data 
analysis was used to average signal values of *GM 18940, 
GM 18942, and GM 18943 fSample ID in HapMap lymphocytes 
from Japanese individuals). 
(TIFF) 

Figure S2 30 different HapMap lymphocytes. The catalog ID 
indicates the database number for HapMap lymphocytes from 
Japanese subjects. 
(DOCX) 

Figure S3 Correlation analysis of SNPs and 6-TGN risk ratios in 

IBD patients. Criteria: r>0.4, p-values<0.01. 

(XLSX) 

Figure S4 SNP genotypes of TPMT. TPMT 1 1 SNP genotypes 

were analyzed for 32 IBD patients. 

(DOCX) 

Figure S5 AZA/5-ASA- induced EAI detected with Express- 
Genotyping analysis. Criteria: AZA/5-ASA-induced EAIa2. 
(XLSX) 
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Figure S6 AZA/5-ASA-induced EAI by AB heterozygous 
genotype analysis of ExpressGenotyping data. This hst shows the 
AZA/5-ASA-induced EAI in more than one (AB genotype in >3 
of 30) HapMap lymphocyte. 

(DOCX) 
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